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Abstract 

The derivatives of a Boolean function are defined up to any order. The 
Taylor and MacLaurin expansions of a Boolean function are thus obtained. 
The last corresponds to the ring sum expansion (RSE) of a Boolean func- 
tion, and is a more compact form than the usual canonical disjunctive 
form. For totalistic functions the RSE allows the saving of a large num- 
ber of Boolean operations. The algorithm has natural applications to the 
simulations of cellular automata using the multi site coding technique. 
Several already published algorithms are analized, and expressions with 
fewer terms are generally found. 

Keywords: Boolean derivatives; Cellular automata; Minimization of 
Boolean functions; Multi-site technique. 



1 Introduction 

This work is based on the concept of Boolean derivatives, introduced in the 
context of cellular automata by G. Vichniac Q. A cellular automaton (CA) is 
defined on a lattice by an interaction range (for instance on a two dimensional 
square lattice, with nearest and next to nearest neighbors interactions), and by 
an updating rule that gives the future value (state) of a lattice variable given its 
present state and the state of its neighbors. The rule is applied in parallel on all 
the sites of the lattice, and can be either deterministic or probabilistic. From a 
computational point of view, the simplest case for the rule is a Boolean deter- 
ministic function, and, if not otherwise specified, we shall refer to this situation 
in the following. Cellular automata are often studied from a numerical point of 
view. Generally large lattices and long time simulations are required, and this 
originates the problem of developing efficient algorithms for the simulations of 
cellular automata on general purpose computers and sometimes on dedicated 
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machines. For the first hardware resources, a technique that allows an efficient 
use of the memory and CPU is the Multi Site Coding (MSC) technique [g |, [§. 
This technique implies that the rule is coded only using bitwise operations. Al- 
though standard bitwise expressions (canonical forms) are easy to generate given 
a Boolean function, the minimization of the number of required operations is 
believed to be a np-complete problem Q . In Section 2 we recall the basic defini- 
tions and in the following section we introduce the Boolean derivatives that will 
lead to the Taylor and MacLaurin series for a Boolean function. They are more 
compact expressions than the standard canonical conjunctive and disjunctive 
forms. In Section 4 this technique is farther developed for the particular case 
of totalistic cellular automata, for which the future state of a cell depends only 
on the total number of neighbors that are in a certain state, regardless of their 
position. The symmetries of the problem allow the saving of a large number of 
Boolean operations. In Section 5 the results are applied to some models that 
appeared in literature. Finally, conclusions are drawn in the last section. 

2 General definitions 

Our fundamental set is B\ = {0, 1}. This is called the Boolean set. Higher 
dimensional Boolean sets are indicated as B n = {0, 1}™. By J- n , m we denote 
the set of functions / : {0, 1}" i— > {0, l} m . T n stands for !F n< x- 

To an element a = (ax, . . . , a n ) € B n corresponds a number a £ [0, 2™): 

n 

a = ^a,-2 1 - 1 ; 

i=l 

and to each number a £ [0, 2") corresponds a n-tuple (ax, . . . , a n ) £ B n : 

<Xi = La/2 ,_1 J mod 2, 

where [a\ stands for the integer part of a. An integer number is thus an ordered 
set of Boolean numbers (bits). In order to emphasize its Boolean structure we 
shall refer to these sets with the name of bitarray, whose dimension is that of the 
space in which it is defined. We introduce a partial ordering between bitarrays 
saying that a < b if, for all i, di < bi, where < 1. We can thus substitute the 
expression a G [0, 2™" 1 ) by a < 2" _1 — 1 or simply a S B n . A Boolean function 
/ is called monotone if a < b implies f(a) < f(b). 

This mapping between numbers and Boolean n-ple corresponds to the repre- 
sentation of integer numbers in computers, the integer division by a power of two 
being equivalent to a right shift (in FORTRAN [a/2 i_1 J = SHIFTR(a, i - 1)), 
and the modulo two operation to take the leftmost (less important) bit. 

Let us introduce the most common Boolean operations. If applied to a 
bitarray, they will act bit by bit. There are 2 2 Boolean functions in T n . With 
n = 1 the most important function is the negation (NOT), indicated by the 
symbol -i, or by a bar over the argument. With n = 2 there are 16 functions. 
The ones that exist on (almost) all computers are the AND, OR and XOR 
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operations. The AND (symbol A) gives one only if both the arguments are one 
(it is equivalent to a multiplication of the bits considered as integer numbers), 
the OR (symbol V) gives one if either one or the other argument is one, while 
the XOR (symbol ffi) corresponds to the sum modulo two. Notice that the 
XOR operation accounts both for the sum and for the subtraction or negation 
(a ffi a = and a ® 1 — ->a) . The AND has higher priority than the OR and 
XOR operations. The negation has the highest priority. The OR and the XOR 
operations are distributive with respect to the AND. The XOR operation can be 
expressed by the NOT, AND and OR: a © 6 = ->a A b V a A ->b. Equivalently the 
OR can be expressed by the AND and XOR operations: aVb = a(BbS>aAb. If 
two Boolean quantities a and b cannot be one at once, both the expressions a®b 
and aVb give the same result. In the following we shall emphasize this condition 
by writing a + b, and indeed in this case the usual sum operation can be used. 
On certain computers (namely on the CRAY), the logical and numerical unities 
can work in parallel. By mixing Boolean and arithmetic operations it is possible 
to speed up the actual calculations M. 

Often in the literature Q the conditional negation is indicated by a b with 
the meaning that a = ->a and a 1 = a. This is equivalent to ->(a © b) or to 
a © b © 1 . In this work we prefer to assign a different meaning to exponentiation, 
more similar to the usual one. We define a = 1 and a 1 = a. With this 
definition the expression a b is equivalent to (a © 1) A b ffi 1 or to a V When 
applied to bitarrays, the exponentiation is performed bit by bit and the results 
are afterwards ANDed together, 



a b 



= /\o»' = /\o i V-6 i . (1) 
i=l i=l 

We note that a b is equal to one if and only if a < b. 

A Boolean function / G T n : x G B n —> f{x) is defined by giving its results 
for all the possible (2™) configurations of its arguments (truth table). It is 
possible to obtain an expression for the / only containing the AND, OR and 
NOT operations from this table. It is sufficient to give / _1 (0) or / _1 (1). The 
two canonical forms are 

n 

f( x ) = V A ^{ x i © a i)\ disjunctive form, and 

n 

f(x)= A \f Xi® en; conjunctive form. 

aef-Ho) i=i 

These canonical forms are the standard starting points for the problem of 
minimizing the number of required operations given a function. It is possible 
to demonstrate || that the minimal expression for a monotone function only 
contains the AND and OR operations; the expression is unique and easy to 
compute. 

As the NOT and the OR can be expressed by the AND and XOR operations, 
any function / can be given in terms of the latter two operations (Ring Sum 
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Expansion, RSE) S. This form is identified by a Boolean vector ff,i £ B„ 

/(*)=© /iAx*. (2) 

Since the number of different Boolean vectors fi and of functions / £ T n is 
equal, the RSE is unique. 

3 Boolean derivatives 

Following Vichniac Q, we define the derivative of a Boolean function / 6 T n 
with respect to its i-th argument a;, as 

df_ 



— f (.X\ j • • • j ^i) • • • j ) © f \X 1 ) ■ ■ ■ ! ? • • • ) ) ■ 



This (first order) derivative expresses the dependence of the function by 
its i-th argument: df/dxi is one if / changes when changing Xi, given the 
configuration xi, . . . , Xi-i, Xi+i, ■ ■ . , x n . If the derivative of / with respect to 
Xi is one regardless of the other arguments, than the rule changes its value 
whenever Xi does. In Ref. a rule that shows this behavior is called a toggle 
rule. 

This definition is consistent with the common expectations: the derivative 
of the identity function is one, and the derivative of a constant (0 or I) is zero. 
Moreover, the derivative is linear with respect to the XOR operations, and it 
follows the standard rule for the derivative of a product, 

ox ox ox 
We can extend the definition to higher order derivatives. For example, a 
second order derivative with respect to Xi and xj is defined as 

d 2 f 

f (x \ , ■ ■ ■ , Xi , . . . , Xj , . . . , Xn) f {x\ : . . . : ~ 'X^,..., Xj , . . . , x n ) 

f (x\ , . . . , Xi , . . . , iXj , . . . , X n ) © f \X\ , . . . , ~~ 'Xj , . . . , — iXj , . . . , x n ) . 



dxidxj 



Note that the definition is consistent with the usual chain rule for derivatives, 

i.e., 



d 2 f d [df 



(3) 



dxdy dy \dx / 

A second order derivative with respect to the same argument is identically 
zero. 

We introduce a more compact definition for the Boolean derivatives. Indi- 
cating with x the bitarray of the arguments (xi, . . . , x n ) and with 5 a (constant) 
bitarray of the same dimension, we define 



d s f(x) = ®f(x®a). 



t<8 
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It is immediate to verify that dsf{x) is equal to the partial derivative in x of / 
with respect to the variables that correspond to nonzero bits in 5. For instance, 
indicating with £W a bitarray in which only the i-th bit is set to one (i.e., 

8f — Si j, where the latter is the usual Kronecker symbol), we have 




We can now state our most important result. For a Boolean function the 
Taylor expansion is always finite. Let us start with a perturbation on only one 
variable. If y = x © from the definition (p|) of the derivative we get 

m=m@d 8m f{x). 

Generalizing 

f(x®8) = ®d a f(x), 

with d Q f{x) = f(x). 

Using our definition of the exponentiation (H), we can substitute the XOR 
over a < 5 with a XOR over a 6 B n . We need a test function that gives one if 
a < 6 and zero otherwise, and from the consideration after Eq. (|l|) this can be 
expressed as a b . Finally we obtain 

f(x®5)= 8 a d a f(x). 

As a noticeable consequence, we can expand a function starting from (the 
bitarray that has all the bits to zero), obtaining the MacLaurin series 

/(*)= * Q A/a, (4) 
aEB„ 

where 

f a = 0«/(O); (5) 

Which is the ring sum expansion of the function /, Eq. (|J). 

Let us explicitly write down the formula (f|) for an elementary CA, whose 
evolution rule depends on the cell itself (y) and on its nearest neighbors (x and 
w). Locally 

y' = f(x,y,w) 

where the prime indicates the future value of the cell. The MacLaurin expansion 
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of / is given by 



y' = /(0,0,0)© 
df 



df 

dx 000 8y 

df 

x Ay A 



z A 



dxdy 
x Ay A z A 



0,0,0 

) x A z A 



df 

dz 

df 



0,0,0 

df 



dxdi 



0,0,0 



0,0,0 



> y A z A 



df 



dydz 



0,0,0 



dxdydz 



The first order derivatives of all the elementary CA can be found in Rcf. §. 
Higher order derivatives can be obtained by using the chain rule (^). Otherwise, 
the array of derivatives fi in zero can be obtained from the truth table f(j) via 
the matrix Mij 

f. M id Af(j); 



where 



mod 2. 



The matrix M. can be recursively generated considering that 

M it o = 1; 

M ,i = (j>0); 

Mij = Mi-ij ®Mi-i d -i {i,j > 0); 



(6) 



To show an application of the MacLaurin expansion, let us examine the 
expression normally used to select between two random bits a and b according 
to a third one (r), 

f(r) = rAaV^rAi (4 operations). 

We only consider the explicit dependence of the function on r. To write down 
the ring sum expansion of f(r) we need 

/(0) = b; 
di/(0) = /(0) ©/(!) = &© a. 



The RSE for f(r) is 

/(r) = 6® r A (a © 6) 



(3 operations). 



We consider it a good result to save one operation out of four in such a widely 
used and (apparently) simple expression. Other examples can be found in sec- 
tion 5. 
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4 Totalistic rules 



The power of the algorithm is particularly evident when applied to totalistic 
CAs. The transition rule for these automata depends on the sum of the cell 
values in the neighborhood, 

n 
i=l 

Any totalistic evolution rule can be written as 

/K...,^)=/(TW)=^r fe Axi n) (7) 

k=0 

where Xk ^ 1S onc 

if T< n ) = k and zero otherwise (totalistic characteristic func- 
tions). Only one term contributes in the sums of equations (Q) so that we can 
use the arithmetic summation. The quantities take the value zero or one and 
define the automaton rule. Probabilistic CAs may be implemented by allow- 
ing the coefficients of equation (Q) to assume the values zero and one with 
probabilities pk (see also the last example of Section 5). 

A totalistic function / is completely symmetrical with respect to its argu- 
ments S . This implies that the derivatives of / of same order are all functionally 
equals. In particular, as the derivatives of the MacLaurin expansion (Q) are cal- 
culated in zero, they are actually equals, and thus can be factorized. This leads 
to 

f( Xl , . . . , Xn ) = / (tW) =/„e/iA Ci n) © h a ©■■■/« a ei n) ; (8) 

where the fi represents the derivative of order i of / in (|5|) , and the t£"^ are 
the homogeneous polynomes of degree i in the variables x\, ...,x n (using the 
AND for the multiplication and the XOR for the sum) 

= x 1 ®x 2 ®---®x n , 
£2™' = xi A x 2 © xi A x 3 © ■ • • © x n -i A x n , 

= xi A x 2 A ■ • ■ A x n . 

The functions ^ satisfy some recurrence relations. The first one is based on 
the idempotent property of the AND operation (a A a = a) and the nullpotent 
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property of the XOR operation (a © a = 0) 



£1 : irreducible; 
£2 : irreducible; 

6 = 6 a a; 

£4 : irreducible; 

6 =€4 a a; 

& = & A 

& = & A £3 = & A & A £1; 
^8 : irreducible; 

The second property is based on the separation of the variables in two groups 
(bisection). Let us call X the group of the variables {x\, . . . ,x n ), with L we 
indicate the left part of X up to some index j, and with R the right part of X 

L = {x\, . . . , Xj) 
R , ■ ■ ■ , x n ). 

We have 

£i(X) = &(L) © 6_i(L) A Zi(R) © &_ 2 (L) A^)©... 
©&(L) A&_i(i2)e&(i2). 

As an example, let us explicitly calculate the & for eight variables. We bisect 
homogeneously the set X = (x\, . . . , x s ) first into L, R, and then into LL, LR, 
RL, RR. We have 





= X\ © x 2 


61 {LR) 


— x 3 © X4 


URL) 


= x 5 © x 6 


Hi{RR) 


— x 7 © x s 




= x\ A x 2 


ULR) 


= x 3 A x 4 


URL) 


= x 5 A x 6 


M RR ) 


— x 7 A x$ 



'UL) 


= Si{LL)®Zi{LR), 




UR) 


= €i{RL) © £i{RR); 




UL) 


= ULL)(Bti(LL)At; 1 (LR)Q 


>ULR), 


UR) 


= URL) © URL) A £i(RR) < 


$URR>: 


' UL) 


= 6(i)A6W, 




UR) 


= 6(^)A6(fl); 




UL) 


= 6(ii)A6(ii?), 






= U r l)aURR); 
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^ 8) = z 1 (L)®Zi(R); 

sf } = 6 (£)©&(£) a &(#)©&(#); 

sj 8) =6A£i; 

d 8) = £4 (L) © 6 (i) a & (#) © 6 (£) a 6 © 6 (i) a 6 (#) © £ 4 (i?) ; 

d 8) =e 4 Aa; 

sf } =£4 a 

U 8) = £4(L)A&(JJ). 

/"q\ 

where £1 = £fe(X). Taking into account the common patterns that appear in 

the expressions of an d Ci 8 "* > we only need 34 operations to build up all the 

As) 

Si 

The extension of the calculations to 9 variables, only adds a small number 
(16) of operations 



A9) 

si 


AS) 
— si 


© £9, 


t(9) 


t(8) 


© C (8) i A19 (2 < i < 


A9) 
S9 


t(8) 
— S8 


A x 9 ; 



even though a careful bisection of the set of variables implies fewer (39) opera- 
tions. 

A kind of normalization condition on the £j is given by 



i=l 



and can be used to save operations in building an expression containing a XOR 
of £j . Another useful relation is 



n — 1 n — 1 

An) 



V(^®^+i) = ®r ) - (10) 



i=l z=l 



We now have to build up the derivatives (in zero) of a totalistic function /, 
Eq. (^). There are only n + 1 independent derivatives /j, (i = 0, . . . , n), as all 
the derivatives of the same order i are equal. We have 



/i = ©M,iA/(T), 



T=0 

where the matrix A4 is defined in Eq. 
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For completeness, we report the expressions for the Xk an( ^ Xk' i 



X9) 



(8) 
Xl 


= 6® 


6 ® 


?5 ©€7, 


(8) 

xF 


= 6© 


6 ® 




(8) 
X3 








(8) 
X\ 




6 ® 




(8) 

xl 


= £ 5 ® 


ft, 




(8) 

Xe 


= £ 6 ® 


fr, 




(8) 
X7 


= 6, 






(8) 

Xs 


= £ 8 ; 






(9) 
Xl 


- Y (8) 

— Xi 






(9) 


(8) 




(9) 


X2 


= X2 




Xs = 


(9) 
Xg 


= 







(11) 



(8) 
X 8 % 



Obviously, the xi 9 "* are only formally similar to the x^ an< l they are calculated 



with nine variables. 

We note that the normalization condition on the 



(n) • 
XL is 



E i 
Xfc 

k=0 



(n) 



(12) 



from which Xo^ can be obtained. 

Putting all the stuff together, we need a maximum of 1024 operations for 
a generic rule with eight arguments, and 2304 operations for a generic rule 
with nine arguments (if all the operations are explicitly developed); 50 (resp. 
57) operations for a generic totalistic rule with eight (resp. nine) arguments 
using directly the of Eq. (||) and 73 (resp. 82) using the hk of Eq. (^). 
These numbers should be compared with the ~ 3000 operations of the standard 
disjunctive form for a function of eight arguments whose truth table is half filled 
with ones and with the ~ 600 operations required for a totalistic rule with nine 
arguments as described in Ref. ||. 

We can see that the Taylor expansion of a Boolean function allows a big 
saving if the function itself depends symmetrically on the variables (i.e., it is 
a totalistic function). Sometimes a function depends in a totalistic way only 
on part of the variables (see e.g., the Life rule in the following section). After 
rearranging the indices, 

f {X\ , . . . , Xi , Xi+i , . . . , X n ) — f {^Cj{x \ , . . . , Xi) , Xi^-i , . . . , X n ) , 

where g(x\, . . . , Xi) is a totalistic function. If this happens, we have 
f{...,x h ...,x k ,...) =/(..., x fc ,..., x jt ...) U,k<i) Vx e B n . 
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From a computational point of view, / depends symmetrically on Xj and Xk 

if 

V Ui © /*©*) = 0, (13) 

where 5 = 8^ © . Symmetries among more than two variables can be 
obtained via the transition property. 

5 Some applications 

In this section we shall apply the algorithm to some problems, chosen among the 
ones that appeared in literature. Some of them were investigated with efficiency 
in mind, so they are supposed to be carefully studied with the aim of reducing 
the number of required operations. 

The first example is the totalistic two-dimensional CA M46789 M. The 
future value d of a cell c is determined by the value most prevalent in its Moore 
neighborhood (nearest and next to nearest neighbors, nine variables), with a 
twist in case of a marginal majority or minority. In terms of Eq. (0) the rule is 
defined as 



rk 



1 if £ = 4,6,7,8,9; 
otherwise. 



The twist in the majority provides a kind of frustration that simulates a mobile 
interface according to the Allen-Calm equation pp| . 

From the general expression (^), we get the simplified expression 

d =&[ie£i(i ©&)]©&, 

for a total of 39 operations. 

The second model is the game of Life fli"|| . This well known game has 
recently shown to be a good tool model for the problem of the self-organizing 
criticality |l2|, [l|. The evolution rule for Life depends separately on the sum 
of the eight nearest and next-to-nearest neighbors (outer Moore neighborhood), 
and on the central cell itself. 

The evolution rule can be expressed saying that a dead (zero) cell can become 
alive (one) only if it is surrounded by three alive cells, while survival only occurs 
if the alive cell is surrounded by two or three alive cells. Developing first the 
rule for the central cell c, we get 

c' = xf©xf Ac, 

(8) 

where c represent the updated value of the central cell, and the Xk are calcu- 
lated on the outer Moore neighborhood. 

The substitution of the expressions for the x^ 8 ' ([ll]) and simplification gives 

d = & A [c © (1 © c) A £i] A (1 © U) . 
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As a © -icj A b = a V b we have 



c' = 6 A(lffi&) A(cV&), 

that implies only 33 operations, to be compared with the ~ 170 ones reported 
in Ref. §. 

We can also apply the method to non totalistic rules. Let us examine the 
Kohring rule |L4| for an FHP |H| gas with obstacles. The collision rules are the 
same of the original FHP model, with a set of four body collisions. Let us label 
the six directions in a counterclockwise way with the indices ranging from one 
to six. The operations on the indices are supposed to be modulo six. All the 
Boolean quantities are actually elements of some array of integer words: we do 
not consider here the spatial indices. If the Boolean variable Xi is one, there is 
an incoming particle on the site from direction i. Collisions occur for 

!(1, 0, 0, 1, 0, 0), two particles collisions; 
(1, 0, 1, 0, 1, 0), three particles collisions; 
(1, 1, 0, 1, 1, 0), four particles collisions. 

The index j = 0, . . . , 5 provides for all the rotational invariant configurations. 

After the application of the collision rule, the variable yi equal to one means 
that there is a particle outgoing from the site with direction z + 3. If the particles 
travel unperturbed, the updating rule is just yi = Xi. On each lattice site there 
is an additional bit, indicated by a, to code the local conservation of angular 
momentum. If a is one and a collision takes place, then all the particles on 
the site turn counterclockwise of 7r/3 (i.e., yi — Xi-i), otherwise the rotation is 
clockwise (yi = Xi + i). The a bit is reversed each time a collision takes place. 
Finally, a bit b = (1) indicates the presence (absence) of an obstacle. The 
meaning of b is reversed from that in Ref. |l4j| for convenience. In the case 
b = no collision occurs, but the velocities of the particles are reversed, i.e., 
yi = Xi+%. For further details about the implementation we refer to Khoring's 
work jb§. 

First we want to obtain the expression for a bit c that indicates the occur- 
rence of a collision. The two cases of zero and six particles can eventually be 
included among the collisions, without having any consequence. There are no 
symmetries among the variables in the truth table of the collisions (see Eq. (|l3|)), 
so it is preferable to divide them into two groups, two and four particles colli- 
sions in one group, and three particle collisions in the other one. The first group 
is characterized by symmetries between Xi and Introducing the auxiliary 

variables Wi — Xi © x i+3 (only three of them are really needed), we get 

c(0,2,4,6) = x^ 3) (wl,w2, W 3); 

where c(i,j, . . . ) indicates the contribution to the collision bit by the i, j, . . . 
particles collisions, and the Xu are the totalistic characteristic functions for 
three arguments. From Eq. (|l2|), (^|) and ( |TT| ) we obtain 

-.c(0,2,4,6) = wlVio2Vi«3. 
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Three particles collisions occur when £3, X5) or {x%, £4, x§) are all zero 
or one, that is 



n C (0,3,6) 



Xi 3) (xi,x 3 ,x 5 ) +x^ 3) (a;i,a:3,a;5)l V I Xi'{x 2 , x A , x 6 ) + X2' ( x 2, xq) 

£,\{xi,Xz,Xc,) ©6(^1^3,^5)] V [£i(x 2 ,X4,,Xg) ©6(^2,334,^6)] ■ 

Using the property (|10|) we get 



.( 3 ). „ „ \ , ,.(3)/ 



->c(0, 3, 6) = (xi © x 3 ) V (x 3 © x 6 ) V (#2 © 14) V {x 4 © «e); 
c=-.[-.c(0,2,4,6) A -.c(0, 3, 6)] . 

This expression for the collision bit is equal to that found in Ref. |Q. 

The expression for the yt can be though as a function of a, b, c. Developing 
the expression we get 

Vi = x l+3 Ob A {x l+3 © ^ © Zi), 

where 

Zi = [Xi © Xi+i © (xi+i © Xi-i) A a] Ac. 
Finally, we notice that 

Zi © z i+3 = [wi © w l+1 © © iOi_i) A a] A c; 

but when c = 1 all the Wi are zero, so Zj = 2:^+3. We need one more operation 
to reverse the angular momentum bit in case of a collision, a' = c © a. Taking 
into account the common patterns in the expression of c,Zi and y,, we only need 
35 operations to update all the six velocities, and 14 arrays (six for the old 
values of the particles, six for the new values, one for the angular momentum 
and one for the collision bits). For comparison, in Ref. [H the algorithm needs 
74 operations and 16 arrays. 

Incidentally, we observe that only six arrays are really needed to store the 
configuration, without the needing of a translational phase. Indeed, the rule 
only implies an (eventual) exchange of the particles among the planes (the RAP1 
machine |l(| is based on this consideration). The translation of the particles 
can be taken into account by a logical shift of the origin of the arrays. The 
procedure is still vectorizable, but the mapping between the lattice and the 
arrays of integer words is indeed more complex, so maybe it is not worth doing 
the efforts unless perhaps for a dedicated hardware. 

The last example involves a probabilistic totalistic CA for the simulation 
of the Ising model ||. The rule depends in a totalistic way on the outer Von 
Neumann neighborhood (the north, east, south and west neighbors). The rule 
can be expressed as 

k=Q 



13 



where the are random bits equal to one with predefined probabilities pk — 
p{ru)- Building the from the & as in Eq. ([ll]), we get the quoted result 
of 22 Boolean operations and four arithmetic summations. Writing down the 
RSE (§), we have 

4 

d = r ©0Si 

i=l 

where the Sj are random bits with probability p(si), obtained as 

si = r ©n, 

52 = r © r 2 , 

53 = r © ri © r 2 © r 3 , 
s 4 = r © r 4 ; 

and considering that p(a © 6) = p(a) +p(6) — 2p{a)p{b). 

Using the approach described above, we need nine operations to build up 
the £i and eight operations for c'. 

6 Conclusions 

In this work we extend and complete the notion of the derivatives of a Boolean 
function already introduced in Ref. Q. We are thus able to derive the Taylor 
and MacLaurin series of a Boolean function. The latter represent the ring sum 
expansion for a Boolean function, which is more compact than the canonical 
conjunctive and disjunctive forms. Moreover, for totalistic functions (i.e., for 
functions completely symmetric in their arguments) very compact expressions 
are found. These ideas have wide applications in the development of faster 
algorithms, in particular for cellular automata simulations, and in the design 
of digital circuitry. As examples of physical applications, we analyze already 
published optimized algorithms, involving both deterministic and stochastic au- 
tomata. We found that our procedure generally leads to more compact expres- 
sions. We think that the Boolean derivative is not limited to the minimization of 
Boolean functions. Work is in progress about the connection between Boolean 
derivatives and the chaotic properties of cellular automata, possibly leading to 
the definition of Lyapunov exponents for discrete systems. 
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